CONVEXITY PRESERVING JUMP DIFFUSION MODELS 
FOR OPTION PRICING 
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Abstract. We investigate which jump-diffusion models are convexity 
preserving. The study of convexity preserving models is motivated by 
monotonicity results for such models in the volatility and in the jump 
parameters. We give a necessary condition for convexity to be preserved 
in several-dimensional jump-diffusion models. This necessary condition 
is then used to show that, within a large class of possible models, the 
only convexity preserving models are the ones with linear coefficients. 



1. Introduction 

A model for a set of stock prices is said to be convexity preserving if the 
price of any convex European claim is convex as a function of the underlying 
stock prices at all times prior to maturity. As is well-known, this property is 
intimately connected to certain monotonicity properties of the option price 
with respect to volatility and other parameters of the model. Generally 
speaking, if the option price is convex at all fixed times, then it is also 
increasing in the volatility. This robustness property motivates the study of 
convexity preserving models in finance. 

Although these issues have been studied quite extensively during the last 
decade, compare [3], 0) El and for the case of one-dimensional dif- 
fusion models, 0, |H and for several-dimensional diffusion models and 
[2] for one-dimensional jump-diffusion models, the general picture for more 
advanced models is not yet fully understood. In 0, a sufficient condition 
for the preservation of convexity in one-dimensional models with jumps is 
provided. That condition, however, is not a necessary condition for preser- 
vation of convexity. The main contribution of the present paper is to give 
a necessary condition for convexity to be preserved in jump-diffusion mod- 
els in arbitrary dimensions. We also use this necessary condition to show 
that, within a large class of possible models, the only higher-dimensional 
convexity preserving models are the ones with linear coefficients. 

To analyze the convexity of an option price we employ the characterization 
of the price as the unique viscosity solution to a parabolic integro-differential 
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equation 

(1) u t = Au + Bu 

with an appropriate terminal condition. In this equation, A is an elliptic 
differential operator associated with the continuous fluctuations of the stock 
price processes, whereas B is an integro-differential operator associated with 
the jumps of the stock price processes, compare Section 01 below. Preserva- 
tion of convexity of the solution to the equation Q is dealt with using the 
notion of locally convexity preserving (LCP) operators. This concept was 
introduced and analyzed in JOJ , and also used in [3] and [Sj • Following these 
references, we show that the condition that A4 = A + B is LCP at all points 
is necessary for convexity to be preserved. We also show that A4 is LCP if 
and only if both A and B are LCP, i.e. 

The model is 1 ( . . . T „„ n ( Both A and B 

.. I J M is LCP \ I 

convexity > =>■ < . , > < are LCP 

I at all points I . , 

preserving I v J I at all points 



compare Theorem 14.21 Thus the characterization of LCP models breaks 
down into two easier problems: (i) to describe which diffusion models are 
LCP, and (ii) to describe which jump structures are LCP. Issue (i) has been 
dealt with in 10 and [3], and (ii) is dealt with in Theorem 15.11 below. 

The present paper is organized as follows. In Section |2] we introduce the 
model and we motivate the study of convexity preserving models by means 
of a monotonicity result. In Section |3] we prove a technical regularity result 
which is used in the sequel. In Section |I] we introduce the LCP-condition, 
and we show that a model is convexity preserving only if both the differential 
operator A and the integro-differential operator B are LCP at all points. In 
Section [21 we investigate which jump structures are LCP. This investigation 
is continued in Section H3 for models with only a finite number of possible 
jump sizes, where we show that, within a large class of possible models, 
all convexity preserving models have linear diffusion coefficients and jump 
structures. 



2. The model and a monotonicity result 

We consider a market consisting of n different stocks, the prices of which 
are modeled by an n-dimensional stochastic process X(t). To specify X, let 
W be an n-dimensional Brownian motion, and let v be a Poisson random 
measure on [0, T] x [0, 1] with intensity measure 

q(dt, dz) = X(t) dtdz, 

where A is a deterministic function. Let X be a jump-diffusion satisfying 
the stochastic differential equation 

dX = (3{X(t-),t)dW + [ <p{X(t-),t,z)v{dt,dz). 

Jo 

Here (3 = (Aj)^=i is an n x n-matrix, <p = (<fti, ■■■,4>n) is an n-dimensional 
vector and v is the compensated jump martingale measure defined by 

v(dt, dz) = (v — q)(dt, dz). 
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Remark In this model, jumps occur according to a Poisson process 



with intensity A(i). Associated with each jump is a label z G [0,1]. The 
interpretation is that a jump of Y at time t with label z results in a jump of 
size 4>i(X(t—),t, z) in the zth coordinate of X. Between jumps, X follows a 
continuous diffusion governed by the diffusion coefficient 0(X(t),t) and the 
drift -\f 1 (f>(X(t),t,z)dz. 

We denote by Wf the space (0, oo) n , and we say that a function g : M™ — * 
M. is of at most polynomial growth if there exist constants m and C such 
that 



for all x G Wf . Given a continuous contract function g : M" — > M of at most 
polynomial growth, the price at time t G [0, To] of an option paying g(X(To)) 
at time T G [0,T], is u(X(t),t). Here the function u : R™ x [0,T ] -> R is 
given by 



where the indices indicate that X{t) = x. Note that the conditions (Al)- 
(A7) specified below and the polynomial growth of g implies that all mo- 
ments of X(T) are finite, compare Section 7 in [Jj. Consequently, the option 
value u is finite. 

Remark We do not address the issue of how to choose an appropriate 
pricing measure, but we rather assume that the model is specified directly 
under the measure used for pricing options. Also note that there is no 
discounting factor in the definition of the option price. Thus we implicitly 
assume, without loss of generality for our purposes, that all prices are quoted 
in terms of some bond price. 

We will throughout this paper work under the regularity and growth 
assumptions (A1)-(A7). When specifying these, D is a positive constant 
and the Holder exponent a is a constant between and 1. 



(Al) For all i,j = 1, ...,n, /% : R+ x [0,T] -> R is in Cl'° (R£ x [0,T]). 
(A2) \eC a ([0,T}). 

(A3) Z]=i\Pij(x,t)\ 2 + \<t>i(x,t,z)\ 2 <Dxf for alH = l,..,n. 

(A4) \P(x,t) - (3(y,t)\ + \<f>(x,t,z) - <f>(y,t,z)\ <D\x-y\. 

(A5) The matrix 0(x,t) is non-singular for all (x,t) G R™ x [0, T]. 

(A6) <j) : R% X [0, T] x [0, 1] -» R n is measurable, and •, z) G C„ ,0 (R r | x 
[0, T]) with the Holder norms being uniformly bounded in the z- 
variable. Moreover, for all x and t we have 4>(x,t,z) ^ for almost 
all z. 

(A7) There exists 7 > — 1 such that 4>i(x,t,z) > 7X4 for i = 1, ...,n. 

Definition 2.1. A model A, (f>) is convexity preserving on [0, T] if for all 

t and Tq with < t < Tq < T , the price u(x, t) of the option with pay-off 
g(X(To)) at Tq is convex in x for any convex contract function g of at most 
polynomial growth. 




g(x)\<C(l + \x\ m ) 



u[x,t) := E X)t g{X(T )) 
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The main reason for studying preservation of convexity is, as mentioned in 
the introduction, that convexity implies certain monotonicity properties of 
the option price with respect to the parameters of the model. The following 
result can be proven in a similar way as Theorem 5.1 in [5], in which the 
one-dimensional case is treated. 

Theorem 2.2. Let two models be given with parameters ((3, A, (f>) and ((3, A, 4>), 
respectively. Assume that 

(i) A(t) < \(t) for allte [0,T]; 

(ii) for each fixed (x,t) € R™ x [0, T] we have cj)(x,z,t) = k(z)(p(x,t, z) 
for some k(z) € [0, 1]; 

(iii) for all (x, t) e R" x [0, T] we have (3(3* < (3(3* as quadratic forms 
(here (3* denotes the transpose of (3). 

Also assume that at least one of the two models is convexity preserving. 
Then, for any convex contract function g of at most polynomial growth we 
have 

u(x, t) < u(x, t) 

for all (x,t) G R™ x [0, To], where u and u are the two option prices corre- 
sponding to the two different models. 

Remark Note that the most important special case of (ii) is when for all 
x, t and z there exists an i G {1, n} such that (f)j(x, t, z) = 4>j{ x i t,z) = 
for all j 7^ i (i.e. the case when at most one component of X and one 
component of X jump at each given time), and 4>i(x,t, z)/<f>i(x, z,t) > 1 if 
4>i(x, z, t) ^ 0. Also note that the condition (iii) is the same as the one used 
for diffusion models in higher dimensions, compare [3]. 

Since all one-dimensional diffusion models and all geometric Brownian 
motions (not necessarily one-dimensional) are known to be convexity pre- 
serving, see , (HI , [H] or |5| , the following consequence of Theorem 12.21 is 
immediate. It is the higher-dimensional analogue to a result in pQ. 

Corollary 2.3. Assume that a model (/3, A, 0) and a convex contract func- 
tion g of at most polynomial growth are given. If n > 2, also assume that 
(3 is the (possibly time- dependent) diffusion matrix of a geometric Brownian 
motion, i.e. (3ij(x,t) = Jij(t)xi for some deterministic functions "fij. Then 
a lower bound for the corresponding option price is given by the option price 
in the model (/?, 0, (ft) with no jumps. 

3. Regularity of the value function 

Under weak conditions, see for example ^2j, the pricing function u is the 
unique viscosity solution of a parabolic integro-differential equation 

(2) u t + Mu = 

with terminal condition 

u(x,T ) = g(x). 
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In this equation, the operator A4 = A+B, where the second order differential 
operator A and the integro-differential operator B are given by 

n 

Au(x,t) := aij(x,t)u XiXj (x,t) 

and 

(3) Bu(x, t) := X(t) / [u(x + 4>(x, t, z), t) — u(x, t) — (ft(x, t, z) ■ Vu(x, t)) dz 
Jo 

respectively, and O'ij £1X6 the coefficients of the matrix (3(3* / 2. Under the 
assumptions (Al)-(A7), u is not merely a viscosity solution to (|2|), but it is 
also a classical solution. Indeed, we show below that the value function u is 
regular. The proof has certain similarities to the proof of Theorem 3.2 in j2j 
in which the one-dimensional case is treated. However, for the convenience 
of the reader, and since the proofs differ at some points, we include the 
higher-dimensional version in its full detail. 

Theorem 3.1. Assume that g £ C^(R") and that g is globally Lipschitz 
continuous. Then u G C^QR^ x [0,T ]). 

Proof. Let ip : 1R+ -*lbea smooth function with ■0' > such that 

s if s > 2 
-1/s if s < 1. 

It follows from Ito's lemma that the n-dimensional stochastic process Y(t), 
where 

Y i (t) = ift(X i (t)), 

satisfies 



000 



n „! 

dY i = b i (Y(t-),t)dt + Y / P ij (Y(t-),t)dW j + \ 

3 = 1 J ° 

on R n x [0,T]. Here 

h(y,t) = U'^rHy^jz^M' 1 ^)^) 



4>i(Y(t-),t,z)v(dt,dz) 



l 



+A(t) J (<k(v,t,z) -i>'{ij- 1 {y i ))ct>i(i'' 1 (y),t,z^ dz, 
P l:j (y,t)=ij'^- 1 (y l ))p ij (i;- 1 (y),t), 



and 



(f>i(y,t,z) = ift (ift 1 (yi) + ((>i(ip 1 (y),t,z)) -y h 
where := •••> i J ~ 1 (yn))- Now it is straightforward to check 

that b, (3 and (ft together with the initial condition g(y) := g(ip~ 1 {y)) satisfy 
the conditions (2.2)-(2.5) in |12j . According to Theorem 3.1 in ^2]) the func- 
tion v(y,t) := u(-0~ 1 (y), t) is a viscosity solution of the integro-differential 
equation 

J v t + E?j=i + Ei=i - A Jo 1 & dz ) % + ^ = 

1 u(y,T ) = g(y), 



(4) 
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where 

h(y,t) = X(t)J (v(y + 4>,t)-v(y,t)^ dz 

and dij are the coefficients of the matrix (3/3* /2. 
Proposition 3.3 in ^21 yields the estimate 

(5) \v(y,t)-v(y,t)\ < C [{I + \y\)\t - i\ 1 ' 2 + \y - y\) 

for some constant C . Together with the assumptions on <p, this implies that 
h G C a (R n x [0,T ]) n C pol (W n x [0,T ]). Consequently, Theorems A.14 
and A. 18 in JU] ensure the existence of a unique classical solution w G 
C po i(R n x [0,T ]jnC , a' 1 (M n x [0,T ]) to equation ©. This classical solution 
io to Q can also be represented (through the Feynman-Kac representation 
theorem) as 

w(y, t) = Ey >t (£° h(Z(s), s) ds + g(Z(T )) 
where Z = {Z\, Z n ) is the continuous diffusion process given by 

dZ i= (bi{Z{t),t) - A(t) MZ{t),t,z)dz^ dt + J2Pij{Z{t),t)dW j 

and Z(t) = y. Since h is Lipschitz continuous in y, it follows from Lemma 3.1 
in ^2] that w is Lipschitz continuous in y, uniformly in t. From the unique- 
ness result Theorem 4.1 in 121 we deduce that v = w. Consequently, 
v G C po i(W l x [0, T ])nCa' 1 (M n x [0, T ]), and therefore h G Cl'°(W n x [0, T ]). 
Applying Theorem A. 18 in ^U] to equation @ once again we find that 
v = w G Cq,' (R n x [0,To]). Transforming back to the original coordinates 
we get u G Ca'^R™ x [0, T ]). □ 

4. The LCP condition as a necessary condition for 

PRESERVATION OF CONVEXITY 

Following ^Jj , see also [1] and |Sj , to investigate which models are convex- 
ity preserving we introduce the notion of locally convexity preserving (LCP) 
models. 

Definition 4.1. An operator D, where T> equals either M, A or B specified 
above, is LCP at a point (x,t) G WL x [0, T] if for any direction v G M. n \ {0} 
we have that 

d 2 v (vf)(x,t)>o 

for all convex functions f G C^(R") n Cp oi (R+ ) with f vv (x) = 0. 
We then have the following key result. 

Theorem 4.2. Consider the following statements: 

(i) The model is convexity preserving. 

(ii) The operator M = A + B is LCP at all points (x, t) G x [0, T] . 

(iii) The operators A and B are both LCP at all points (x, t) G WL x [0, T] . 

We have that (i) =^> (ii) <^=> (iii). 
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Remark Under a few additional growth conditions on (3, 4> an d their deriva- 
tives, it is possible to prove also (ii) ==> (i). Thus all three statements in 
Theorem 14.21 are equivalent, see Theorem 4.3 in [3] for the one-dimensional 
case. We do not pursue this further since we only use LCP as a necessary- 
condition in the analysis below. 

Proof. To prove (i) (ii) we argue as in the proof of Lemma 3.3 in 
Choose (x , T ) G M™ x [0, T] and let g G ^(M^)nC^(]R™ ) with g vv (x ) = 
for some direction v. Define u to be the solution to 

u t + Mu = 

on M™ x [0,T ] with terminal condition u(x,T ) = g(x). Let g G C*(R!J.) be 
convex, Lipschitz continuous and satisfy g = g inside a box which contains 
xq and all possible values of xq + <f>. It follows that Bg = Bg in a (spatial) 
neighborhood of x, so d%(Bg)(xo,To) = d%(Bg)(xo,To). Now, let u be the 
solution to 

u t + Mu = 

on x [0, To] with terminal condition u(x,Tq) = g(x). Then, since the 
model is convexity preserving, u is convex in x at all times prior to To, and 
in particular u vv (xo,t) > 0. Consequently, dtU vv (xo,To) < 0, so it follows 
from Theorem 13 . 1 1 that 

< -d t u vv (xo,T ) = -d^u t (x ,T ) = d%Mu(x ,T ) 
= d 2 v (Mg)(x ,T ) = d 2 v (Mg)(x ,T ), 

which finishes the proof of (i) ==>■ (ii). 

The implication (iii) ==>■ (ii) is immediate from the definition of LCP. It 
remains to show that (ii) ==>■ (iii), i.e. that if M. is LCP, then both A an B 
are LCP. This follows from Lemmas 14.31 and 14.41 below. □ 

Lemma 4.3. Let t G [0,T] be given, and let f G C*(R%) n C^QRIJ.) be a 
convex function with f vv (xo) = at some point xq and for some direction v. 
Then, for any e > there exists a convex function h G C^(M") PI Cp oi (IR™) 
with h vv (xo) = such that d%(Ah)(xo,to) = and \d%(B(f — h))(xo,to)\ < 
e. Consequently, if M is LCP, then also B is LCP. 

Proof. Let (x ,t ) G x [0,T] be given, and assume that / G C^(M^) n 
Cp oi (lR™ ) is a convex function with f vv (xo) = for some direction v. Without 
loss of generality we assume that f(xo) = and V/(xo) = (this can be 
done since Af = Bf = for all affine functions /). Let D be a constant such 
that for all unit directions w we have f ww (x) < D for all x in a neighborhood 
of xq. Let (p : — > [0, 1] be a smooth function such that 

, v f if \y\ < 1 
= { 1 if |y| > 2, 

and let Ci be a constant such that |</?™(y)| < C\ and ^^(y)] < C\ for all 
y G and all unit directions w. Further, let i/j : [0, oo) — > [0, oo) be a 
smooth and non-decreasing function satisfying 

" for s G [0,1/2] 
V>(«) = { Ms for s G [1,2] 
3M for s > 3, 
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where M is a constant satisfying M > 8C±D. Now, for 5 > 0, define the 
function h = h s by 

r\x-x \/S 

h(x) = f [<p (<5 _1 (:r - x )) (x - xq) + xo) + S 2 ip(s) ds. 

Jo 

Then h is the wanted function for some 6 small enough. Indeed, first note 
that h is if \x — xq\ < 5/2. Consequently, d 2 (Ah)(xo, to) = 0. Moreover, h 
is convex if 5 is small enough. To see that d 2 (B(f — h)) can be made small, 
note that 

d 2 v (Bf)(x ,t ) = A f 1 (d 2 v (f(x + 0)) - fw(xo) - dl{cj> ■ Vf(x ))) dz 
Jo 

= A / (d 2 v (f(xo + 0)) - <f>w ■ V/(x )) dz 
Jo 

where we have used f vv (xo) = and f vw (xo) = for any direction w (the 
latter statement follows from f vv = and the convexity of /). Similarly, 

d 2 (Bh)(x ,t ) = A / (d 2 (h(x + 4>))-4>w Vh(x )) dz. 
Jo 

Thus, since V/(xo) = Vh(xo), 

d 2 (B(f - h))(x ,to) = A ^ d 2 (/(x + 0) - h(x + 4>)) dz 

Jo 

= x d 2 (f{x + (j)) - h(x + 4>)) l {z .\<p\ <3S} dz 
Jo 

/•i / M\/8 \ 

:\4>\>3S} dz 

= h + h. 

Here I\ converges to as 5 goes to since (f)(x, t, z) / for almost all z 
by (A6) and d 2 {h(x + 0))1{ Z: |^|<3<5} is bounded uniformly in 5 at x = xq. 
Similarly, 

h = A f 1 {Sm^M/S) + 1^1^(1^1/5)) 1 { ,:| |>3«5} dz 
Jo 

= A / 5\(p\ vv 3Ml {z . l<pl > 3S} dz, 
Jo 

so it follows from (A6) that also I2 converges to 0. Thus, by choosing S small 
enough, we find that d 2 (Bh)(xo, to) is arbitrarily close to d 2 (Bf)(xo, to). □ 

Lemma 4.4. Let t G [0,7] be given, and let f G C^(M^) n C 2 ol (Rl) be 
a convex function with f vv (xo) = at some point xq and for some direc- 
tion v. Then, for any e > there exists a convex function h £ C^(1R") n 
C 2 ol (Wl) with h vv (x ) = such that d 2 (Ah)(x ,to) = d 2 (Af)(x ,to) and 
\d 2 (Bh)(xo,to)\ < e. Consequently, if M is LCP, then also A is LCP. 
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Proof. Without loss of generality, assume that /(xq) = and V/(xo) = 0. 
As is seen in the proof of Theorem 15. II below. 

(6) dl(Bh){x G ,t Q ) = J (h m {xQ + 4>) + 2(j) v -Vh v {x G + (j)) 

+4> v Hh{xo + 



+<j) m ■ (Vh(x +4>)- Vh(xo))J dz 

provided h is convex and h vv (xo) = (here Hh denotes the Hessian of h). 
Thus it suffices to find 5 > and h satisfying h = f on {|x — xq\ < 5} and 
such that Vh and Hh are small on {|x — xo| > 5}. 

To do this, let C\ be a constant such that f ww < C\ in a neighborhood of 
xq and for all unit directions w. Let (p : W 1 — > R be a smooth non-negative 
function satisfying < (p < 1 such that 



ip(x) 



1 if \x — xo\ < 1 
if|x-x |>2. 



Further, let C2 be a constant such that < C2 and |^to| < C2 for all 
unit directions w. Let M > 6C1C2 and let ^ : [0, 00) — > [0, cxo) be a smooth 
and non-decreasing function satisfying 

for s € [0,1/2) 
V>(s) = <( Ms for s G (1,2) 
for s £ (3, 00). 




Now let 



h\x) := / (xV((x - x )/<5) + 5 2 I ' ° f rp(s) ds. 







Then h = h s is the wanted function for some 5 small enough. Indeed, first 
note that Af(xo,to) = Ah(xo,to) since h = f for \x — xo| < 5. Also note 
that straightforward calculations yield that h is convex if 5 is small enough. 
Moreover h{x) = k\x — xq\ + b for \x — xo\ > 35, where k = 3M5 and b is 
some constant. If it? is a unit vector with the same direction as </>, i.e. if 
4> = \<fi\w, then it follows that 

rW 

h Xl (x + 4>) - h Xi (x )= / K iW (x + sw) ds 

Jo 

can be made arbitrarily small (when varying 5) since h XiW is bounded inside 
|x — xo I < 3(5 (uniformly in 8) and vanishes outside this region. Thus the 
last term in the right hand side of (JSJ can be made arbitrarily small. 

Moreover, examining the first three terms of © one finds that these 
together form a second derivative of h, evaluated at xo + (ft, in the direction 
w := v + <p v . Now 

< h ww {x Q + 4>) < k\w\ 2 /\(t)\ 

for > 38. Since k is linear in 5 it follows that also the three first terms 
in © can be made arbitrarily small when decreasing 5. □ 
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5. A CHARACTERIZATION OF LCP MODELS 

In [I] it is shown that, within a large class of models, the only differential 
operators of the form A which are LCP in dimension n > 2 are the ones 
corresponding to geometric Brownian motions. In that sense, there are not 
very many convexity preserving diffusion models in higher dimensions. In 
this section we study the LCP-condition for the operator B corresponding 
to the jump part of X. 

Let H f denote the Hessian of a function /. The following theorem gives 
a precise description of which jump structures eft give rise to an integro- 
differential operator B which is LCP. 

Theorem 5.1. The operator B is LCP at a point (x, t) if and only if for all 
directions v G M n \ {0} we have 

(7) J [f m {x + 0) + 2<f> v ■ Vf v {x + <f>) + 4> v Hf(x + 4>)4>% 

• (V/(x + (/>)- V/(x))) dz > 

for all convex functions f G C*(R!J.) n C^(R™ ) with f vv (x) = 0. 

Proof Assume that /eC*(l n + )n Cp «( M +) is convex and that f vv = at 
x. Straightforward calculations yield that 

d 2 v (Bf)(x) = J 1 ^f vv ( x + £) + 2<j )v -Vf v (x + <j ) ) + <p v Hf(x + <j ) )<Pl 

+4>w ■ V/(x + 4>) - (j) vv ■ Vf(x) - 2<j) v ■ V/o - (f) ■ ^fw^j dz. 

From the assumption f vv = at x it follows, due to the convexity of /, that 
V/u = V fw = at x. Thus 

d 2 v (Bf)(x) = (f m {x + </>) + 24> v -Vfvix + ^ + ^Hfix + ^l 

+(/> vv -{Vf(x + <P)-X7f(x))^dz. 

Consequently, B is LCP if and only if (JJJ) holds for all directions v € W 1 \ {0} 
and all / G C*(Wl) H C 2 pol (Rl) with f vv (x) = 0. □ 

Corollary 5.2. Let (x,t) G R" x [0,T\. If 

(8) f <p vv ■ (Vf(x + <p)- V/(x)) dz > 
Jo 

for all convex functions f G C^(IR™) n C 2 ol (W^_) and all directions v G R n \ 
{0}, then the operator B is LCP at (x,t). 

Proof. Since / is convex, and since f vv + 2(j) v ■ S7f v + 4> v Hf(j)* is the second 
derivative of / in the (v + <^)-direction, it is clear that (JHJ) is sufficient for 
the LCP-condition. □ 

Remark If for all i = 1, ...,n, the function 0, is convex in x at all points 
(x, t, z) where <f>i(x, t, z) is positive, and <\>i is concave in x at all points (x, t, z) 
where 4>i(x,t,z) is negative, then © is clearly satisfied. This sufficient 
condition for preservation of convexity was used in 5 in a one-dimensional 
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setting. Also note that it is possible to show that the condition (JHJ) is strictly 
weaker than the condition (J7J). 

6. The case of only finitely many possible jump sizes 

In this section we investigate models with only finitely many possible 
jump sizes at each time. More specifically, we assume that for each fixed x 
and t, the function z *— » (ft(x, t, z) takes at most finitely many values. 

Theorem 6.1. Assume there are only finitely many jump sizes, and let 
(x, t) G Wl x [0, T] . Then the following conditions are equivalent: 

(i) B is TCP at (x,t). 

(ii) (0) holds for all directions v and all convex f G C*(H.") n C^ ol (W\. ) 
with f vv = 0. 

(iii) |3j) holds for all directions v and all convex f G C^(R") n C^iW^ ). 

Proof. In view of Theorem 15.11 and Corollary 15.21 we only need to show the 
implication (ii) ==? (iii). 

To do this, let / G C*(R%) n C^(M") be convex. Since there are only 
finitely many possible values of (ft, we can deform / to be flat around all 
possible values of x + (ft and also in a neighborhood of x without altering 
the first derivative at these points. Accordingly, the first three terms of the 
integrand in vanish, whereas the last term remains unchanged, which 
finishes the proof. □ 

Theorem 6.2. Assume there are only finitely many jump sizes and that 
the model is convexity preserving. Then, for each i = l,...,n, the "jump 
volatility" 



of the ith asset is increasing as a function of Xi at each fixed time t. 

Proof. Since the model is convexity preserving, it follows from Theorem 14.21 
and Theorem l6.1l that the inequality © holds at all points and for all convex 
functions /eC^)n C% ol (R1). Choosing / = x\ in © gives that for all 
i = 1, n and for all directions v we have 




(9) 




at all points (x,t). Fix an i and let 



ip(x,t) : 




Differentiating ip with respect to Xi gives 
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Using Q with v = e Xi (here e Xi denotes the ith unit coordinate vector), 
integration by parts, 4>i - ► as x% — > and Jensen's inequality we find 

(10) < x\ \ \ (cf)i) XiXi 4>idzdxi 

Jo Jo 

xl4>i{(j)i) Xi -xj I (cj)i)l.dxi ) dz 



< 



(xj4>i(4>i) Xi - Xi$) dz. 



This shows that ip Xi is non-negative, which finishes the proof. □ 

To the best of our knowledge, not very many models in finance have 
increasing volatilities. Instead, models have typically large volatilities for 
small values of the underlyings. If we restrict our attention to these typical 
models, we show below that preservation of convexity is a rather special 
property. 

Theorem 6.3. Let n > 2. Assume that there are only finitely many possible 
jump sizes and that for alii and j, (3ij is a function merely of Xi andt. Also 
assume that for all i = 1, n and for each fixed time t the "total volatility" 



l (5Uxi,t) + ... + Pl(xj,t) + A 2 (t) /J tfjx, t, z) dz 



of the ith asset is not an increasing function of Xi, unless it is constant. If 
the model is convexity preserving, then (3ij and are linear in Xi for all i 
and j, and 4>i does not depend on Xj for j ^ i. More explicitly, there exist 
functions 7^ : [0, T] — ► R and 7^ : [0, T] x [0, 1] such that 

(11) Pij(xi,t) = xaij(t) 

and 

(12) <f>i(x,t,z) = xai(t,z) 
for almost all z. 

Proof. First note that according to Theorem 14.21 both the operators A and 
B, corresponding to the diffusion part and the jump part of X, respectively, 
are LCP. Now note that if the total volatility is strictly decreasing in some 
interval, then either a "diffusion volatility" 



or a "jump volatility" 



l pf l (x u t) + ... + pl(x u t) 

X 1 ; 



t \2(t)f Q 1 tf(x,t,z)dz 



is strictly decreasing in some interval. However, the proof of Theorem 2.3 
in [I] implies that all diffusion volatilities are increasing, and Theorem 16.21 
above implies that all jump volatilities are increasing. Consequently, all 
diffusion volatilities and all jump volatilities are constant in x^. 
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It then follows from the proof of Theorem 2.3 in [|| that all fyj are linear 
in Xj. Moreover, if the ith jump volatility is constant in Xj, then the corre- 
sponding inequalities in lfTU|) reduce to equalities. Since Jensen's inequality 
reduces to an equality if and only if the integrand is constant, we find that 
for almost all z the function fa has to be linear in X{. It thus only remains 
to show that fa does not depend on Xj, j ^ i. To do this we fix j ^ i, and 
we plug v = e x . + se Xj and fa = ji(x,t, z)xi into the inequality ©, where 
s G E and x = {x\, Xj+i, x n ). We find that 

0< / fa{fa) m dz = 2sxi / ji( r Yi) Xj dz + s 2 Xi 7i(ji) XjX:j dz. 
Jo Jo Jo 

Since this expression is non-negative for any choice of s € K, we must have 

(13) f ji(-fi) Xj dz = 0. 

Jo 

Now let x'j G M+. Performing similar calculations as in (Tml) . but with 
v = e r . , we find that 



< (xj-x'j) 2 j (fa(x,t,z)(fa) Xj (x,t,z) - fa(x,t,z)(fa) Xj (x',t,zj) dz 
(fa(x,t,z) - fa(x ,t,z)) 2 dz, 



\Xj Xj\ 



where x' = (x%, Xj-i, x'j, Xj+i, x n ). In view of (fT3)) . the first integral 
vanishes. Consequently, fa is z-almost surely constant in Xj, which finishes 
the proof. □ 

Remark Note that the models satisfying (jlljl and (j!2j) are all convexity 
preserving. Indeed, by explicit solution formulas, 

{ n To n Tq 1 

where 



Ji(T ) =exp|- ( ° ( \{t) ll (t,z)dzdt\ J] II (1+7*W) ? 

L J0 J0 J o<z<10<t<T 



■)v(z,t) 

- 0<i<T 

so Xj(To) is linear in the starting value Xj(0). Consequently, 

u(x,0) = E Xt0 g(X(T )) 
is convex in X(0) provided g is convex. 

References 

[1] N. Bellamy and M. Jeanblanc. Incompleteness of markets driven by a mixed dif- 
fusion, Finance and Stock. 4, 209-222 (2000). 

[2] J. Bergenthum and L. Ruschendorf. Comparison of option prices in semimartingale 
models. Preprint (2004). 

[3] Y.Z. Bergman, B.D. Grundy and Z. Wiener, General properties of option prices, 
J. Finance 51 (1996), 1573-1610. 

[4] E. Ekstrom, S. Janson and J. Tysk, Superreplication of options on several under- 
lying assets, J. Appl. Probab. 42 (2005), 27-38. 



14 



ERIK EKSTROM AND JOHAN TYSK 



[5] E. Ekstrom and J. Tysk, Properties of option prices in models with jumps, 
arXiv.org math. AP/0509232 (2005), submitted for publication. 

[6] N. El Karoui, M. Jeanblanc-Picque and S. Shreve, Robustness of the Black and 
Scholes formula, Math. Finance 8 (1998), 93-126. 

[7] I. Gihman and A. Skorohod, Stochastic Differential Equations (1972), Ergebnisse 
der Mathematik und ihrer Grenzgebiete, Band 72. Springer- Verlag, New York- 
Heidelberg. 

[8] D. Hobson, Volatility misspecification, option pricing and superreplication via 

coupling, Ann. Appl. Probab. 8 (1998), 193-205. 
[9] S. Janson and J. Tysk, Volatility time and properties of options, Ann. Appl. 
Probab. 13 (2003), 890-913. 
[10] S. Janson and J. Tysk, Preservation of convexity of solutions to parabolic equa- 
tions, J. Diff. Eqs. 206 (2004) 182-226. 
[11] M. Kijima, Monotonicity and convexity of option prices revisited, Math. Finance 
12 (2002), 411-425. 

[12] H. Pham, Optimal stopping of controlled jump diffusion processes: a viscosity 

solution approach, J. Math. Systems Estirn. Control 8 (1998) 1-27. 
E-mail address: ekstrom@maths.manchester.ac.uk, Johan.Tysk@math.uu.se 



